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>  _  \  ABSTRACT  C  .....  *- 

\ 

The  Cauchy  problem  for  the  Korteweg-deVries  equation  (KdV  for  short) 

C  qt(x,t)  +  qxxx(x,t)  -  6q(x,t)qx(x,t)  =  0 
(*)  / 

1  q(x,o)  «  q(x )“- 

■is  solved  classically  under  the  single  assumption 

r  n  +  | X I  )4|Q(x)|dx  <  00 

_0O 

-for  t  >  0  via  the  so-called  "inverse  scattering  method".  This  approach, 
originating  with  Gardner,  Greene,  Kruskal,  and  Miura  [9],  relates  the  KdV 
equation  to  the  one-dimensional  Schrodinger  equation: 

(**)  -f " (x,k)  +  u(x)f (x,k)  =  k2f(x,k). 


1  By  considering  the  effect  on  the  scattering  data  associated  to  the 
Schrodinger  equation  (**)  when  the  potential  u(x)  evolves  in  t  according 
to  the  KdV  equation  (*),  one  obtains  a  linear  evolution  equation  for  the 
scattering  data.  The  inverse  scattering  method  of  solving  (*)  consists  of 
calculating  the  scattering  data  for  the  initial  value  Q(x),  letting  it  evolve 
to  time  t,  and  then  recovering  q{x,t)  from  the  evolved  scattering  data. 

Recently,  P.  Dei ft  and  E.  Trubowitz  [7]  presented  a  new  method  for 
solving  the  inverse  scattering  problem  (obtaining  the  potential  from  its 
scattering  data).  Our  solution  of  the  KdV  initial  value  problem  uses  this 
approach  to  construct  a  classical  solution, under  the  assumption  stated  above. 
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SIGNIFICANCE  AND  EXPLANATION 


The  Korteweg-deVries  equation  (KdV  for  short)  arises  as  an  approximation 
in  many  problems  involving  non-linear  dispersive  waves  and  has  been  exten¬ 
sively  studied  in  recent  years  [ 1 , 3-6, 9, 1 1-1 3, 1 5-1 9, 21-28] .  One  approach  to 
solving  this  equation  in  a  more  or  less  explicit  fashion  is  the  inverse 
scattering  method  of  Gardner,  Greene,  Kruskal,  and  Miura  [9],  which  relates 
the  KdV  equation  to  a  one-parameter  family  of  one-dimensional  Schrodinger 
operators.  We  solve  the  Cauchy  problem  for  the  KdV  equation  by  this  method, 
using  the  inverse  scattering  theory  of  Deift  and  Trubowitz  [7],  Previous 
authors  [5,6,9,22]  have  done  this  using  Faddeev's  version  of  inverse 
scattering  [8] ;  our  approach  constructs  a  classical  solution  under  less 
restrictive  conditions  on  the  initial  data.  In  particular,  no  smoothness  is 
assumed.  Some  aspects  of  the  asymptotic  behavior  of  the  solution  are  also 
discussed. 
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CLASSICAL  SOLUTIONS  OF  THE  KORTEWEG-deVRIES  EQUATION 


FOR  NON-SMOOTH  INITIAL  DATA  VIA  INVERSE  SCATTERING 
Robert  L.  Sachs 

1 .  Introduction. 

In  this  paper,  we  construct  classical  solutions  of  the  Cauchy  problem  for 
the  Korteweg-deVries  (KdV  for  short)  equation: 

r  qt  +  "xxx"  *NX  =  °'  >  °'  x  e  K 

(*)  / 

q(x,0)  =  e(x) 

via  the  inverse  scattering  method  of  Gardner,  Greene,  Kruskal,  Miura,  and 
Zabusky  (see  [9,17]).  This  approach  relies  on  the  intimate  connection  between 
the  KdV  equation  and  the  one-dimensional  Schrodinger  equation  with  potential 
u(x) ; 

(**)  -f " (x,k)  +  u(x)f(x,k)  =  k2f(x,k). 

More  precisely,  if  one  considers  the  one  parameter  family  of  Schrodinger 
operators  (parameterized  by  t)  whose  potentials  satisfy  the  KdV  equation, 
then  these  operators  are  unitarily  equivalent  to  one  another.  Moreover,  by 
considering  the  so-called  scattering  data  associated  to  these  operators,  one 
obtains  a  system  of  uncoupled,  linear  equations.  The  inverse  scattering 
method  of  solving  the  KdV  equation  consists  of  first  calculating  initial 
values  for  the  scattering  data  from  the  initial  potential  Q(x),  then  solving 
the  linear  equations  for  the  evolution  of  the  scattering  data,  and  finally, 
recovering  the  potential  q(x,t)  at  time  t  from  its  scattering  data. 
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The  first  two  steps  in  this  procedure  are  relatively  straightforward; 
most  of  the  technical  difficulty  in  applying  the  method  arises  in  the  process 
of  recovering  the  potential  at  time  t.  Previous  authors,  including  Tanaka 
[22]  and  Cohen-Murray  [5],  used  t  a  inverse  scattering  theory  of  Faddeev  [8] 
to  complete  this  step  under  suitable  hypotheses  on  the  initial  data  (involving 
decay  of  several  derivatives).  In  a  later  paper,  Cohen-Murray  [6]  considered 
initial  data  which  was  continuous  with  a  piecewise  continuous  first  deriv¬ 
ative.  Assuming  a  certain  decay  rate,  she  proved  the  existence  of  a  smooth 
solution.  Our  result  does  not  assume  a  specific  decay  rate  but  rather 
requires  that  the  initial  data  belong  to  a  weighted  L1  space.  We  make  no 
assumptions  on  derivatives  of  the  initial  data.  In  recovering  the  potential 
from  its  scattering  data,  we  use  the  recent  work  of  P.  Deift  and  E.  Trubowitz 
[7]  on  the  inverse  scattering  problem.  One  difference  between  this  theory  and 
Faddeev's  is  that  the  so-called  trace  formula  of  [7]  (see  also  Newell  [19]) 
expresses  the  potential  itself  directly  in  terms  of  the  scattering  data.  In 
this  way,  the  solution  q(x,t)  of  the  KdV  equation  (*)  and  the  solution 
h(x,t)  of  the  linear  problem: 

(1.1) 

^  h(x,0)  =  H(x ) 

are  compared  directly,  where  H(x)  is  related  to  the  scattering  data  of 

Q(x).  Our  analysis  consists  in  large  part  of  solving  (1.1)  in  various 
1  " 

weighted,  local  L  and  L  spaces  and  then  using  the  link  between  h(x,t) 
and  q(x,t)  to  extend  these  properties  to  q(x,t). 

Several  other  approaches  to  the  KdV  initial  value  problem  have  been 
developed  besides  the  inverse  scattering  method  used  here.  Saut  and  Temam 
[21]  used  a  parabolic  regularization  of  the  KdV  equation  to  establish  the 
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existence  and  uniqueness  of  solutions  in  the  Hilbert  spaces  Hs  for 
s  >  2.  Bona  and  Smith  [4]  obtained  existence,  uniqueness,  and  continuous 
dependence  on  initial  data  in  Hs  for  all  integers  s  >  2  using  a 
regularizing  term  of  order  3  (the  so-called  regularized  long  wavelength 
equation).  Bona  and  Scott  (3],  using  non-linear  interpolation  theory, 
extended  this  result  to  non-integer  values  of  s.  These  results  show  that  in 

the  Hilbert  spaces  Hs,  the  KdV  evolution  is  not  smoothing  in  the  strong 

2 

sense  but  preserves  the  order  of  L  -differentiability.  Recently,  Kato 
[12,13],  using  the  abstract  theory  of  quasi-linear  evolution  equations,  proved 
the  existence  and  weakly  continuous  dependence  on  initial  data  in  Hs,  s  = 


0,1,  and  the  weighted  Sobolev  spaces 


H2r,r  i  H2r  n  L2( ( 1+x2)rdx) 


and 

„s  2  ,  2  _  2 ,  2bx ,  . 

H  n  L^,  where  «  L  (e  dx) 


for  b  >  0,  s  >  0. 

Kato  also  shows  that  in  fact,  for  s  >  3/2,  an  Hs  solution  belongs  to 
2  s+1 

1-  ( ( 0,T] ; ) ,  which  is  a  kind  of  smoothing  effect.  Other  papers  on  the 
Cauchy  problem  for  the  KdV  equation  on  the  line  include  [11,18,15-27].  For  a 
nice  survey  article,  see  Miura  [17]. 

Our  result,  as  well  as  that  of  Cohen-Murray  [5,6],  exhibits  an  interplay 
between  smoothness  and  decay  in  the  pointwise  and  weighted  L1  senses.  We 
remark  that  these  contrasts  arise  in  the  linear  evolution  equation  (1.1)  as 
well  as  for  the  KdV  equation. 
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(1.2) 


Before  describing  our  results,  we  introduce  some  notation: 
^  L^(a)  :=  LP([a,«)) 

<  =  £u(x)  :  <1+|x|)ru(x)  6  LP((a,«))} 

XP  :=  n  LP(a) 
a  finite 


LP  will  denote  LP(  Note  that  XP  is  a  Frl&chet  space  with  the 

r  r  r 

obvious  topology. 

Our  principal  results  are  as  follows: 

Suppose  the  initial  value  Q(*)  G  Then  there  exists  a  unique 

solution  q(x,t)  of  (*)  with  the  properties: 

1  00 

(i)  q(x, t )  SX  ...  .  n  X  for  every  t  >  0  and  5  >  0 

a-3/4-o  a 

and  t  -*•  q(»,t)  is  continuous  in  this  topology  for  t  >  0. 


(ii)  3^3®qC,t)  €  X1 


(with  continuity  in  t) 

(iii)  q(x,t)  ♦  Q(x)  as  t  i  0 

in  i  <5  for  every  6  >  0. 

As  yet,  we  do  not  know  how  to  prove  that  the  solution  constructed  has  the  same 
asymptotic  behavior  as  x  ♦  -•  or  as  t  *  +“  (modulo  soli tons)  as  the 
solution  of  the  linear  problem  (1.1),  although  we  believe  this  to  be  true 
(See  [1)). 

The  basic  results  of  inverse  scattering  theory,  including  sketches  of  the 
approaches  of  Faddeev  (8)  and  Deift  and  Trubowitz  (7),  appear  in  Section  2 
below.  In  Section  3,  the  link  between  the  KdV  equation  and  the  Schrodinger 
equation  is  described  and  the  evolution  of  the  scattering  data  given.  Using 
this  evolution,  the  reconstruction  of  the  potential  by  the  Deif t-Trubowitz 
method  is  discussed  in  Section  4.  Additional  smoothness  and  decay  properties 


-4- 


of  the  proposed  solution  are  developed  in  Section  5,  while  Section  6  contains 
the  proof  that  our  solution  indeed  satisfies  the  KdV  equation. 

It  is  a  pleasure  to  acknowledge  helpful  conversations  with  Percy  Deift 
and  Jerry  Bona  during  the  course  of  this  work  and  the  support  of  an  A.M.S 
Postdoctoral  Research  Fellowship. 


2.  Results  From  Inverse  Scattering  Theory. 

In  this  section,  we  sketch  the  results  of  the  theory  of  inverse 
scatterinq  for  the  one -dimensional  Schrodinqer  equation  [7,8]  to  be  used  in 
studyinq  the  KdV  initial  value  problem.  Where  feasible,  the  proofs  of  these 
results  are  described  briefly.  Our  discussion  is  divided  into  three  parts: 
the  definition  of  scatterinq  data  for  the  Schrodinger  equation;  Faddeev's 
inverse  method  [8];  and  the  Deif t-Trubowitz  approach  [7].  Our  discussion  is 
for  the  reader's  benefit  and  contains  no  new  results. 


A.  The  Forward  Problem:  Defining  the  Scattering  Data;  Basic  Properties. 

Consider  the  Schrodinger  equation 
(**)  -f  "  ( x ,  k )  +  u(x)f(x,k)  =  k2f(x,k) 

on  the  real  line  <  x  <  Jf°  for  k  real  and  non-zero.  Assuming  that  the 


potential  u(x)  is  real-valued  and  satisfies  the  condition: 

(2.1)  Hun  =  r  (1  +  I x |  )  |u(x)  |dx  <“ 

L1  - 

we  construct  the  Jost  solutions  f^(x,k)  of  (**)  as  follows: 


Solve  the  Volterra  integral  equations: 


r  L2ik(y-x)_1 
( x ,  k )  =  1  +  /  - — - 


u(y)m+(y,k)dy 


(2.2) 


-2ik(y-x )  , 
,  , .  .  rx  e  -1 

(x,k)  =  1  +  J  - ttt; - 


u(y)m  (v,k)dy 


for  m^(x,k)  and  then  define 
(2.3)  f±(x,k)  =  e±ikxm±(x,k). 

Condition  (2.1)  arises  when  estimatinq  the  terms  in  parentheses  by  | y-x | , 
obtaining  a  bound  independent  of  k.  Then  f^(x,k)  satisfy  the 
Schrodinger  equation  (**)  and  have  the  asymptotic  behavior 
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(2.4) 


f+(x,k)  ~  « 


ikx 


as  x  +  ¥° 


,  ,  .  .  -ikx  ^  _ 

f  (x,k)  ~  e  as  x  +  -fl0. 


Since  u(x)  anti  k  are  real,  f+(x,k)  =  f^(x,-k)  are  also  solutions  of 
(**).  (2.4)  implies  f+(x,k)  and  f+(x,-k)  are  linearly  independent  for 

k  *  0  (and  similarly  for  f_)  since 

[f +  (x,k ) , f +(x, -k ) )  =  ff _(x,-k) , f_(x,k) )  =  2ik. 

Thus  there  are  functions  T^(k),R^(k)  defined  for  real,  non-zero  k  by  the 


relations: 


(2.5)  i 


1  Vk) 

f  <x,k)  =  =r-7rr  f  ,  (x,-k)  +  — rCT  f  .  (x,k) 


T+(k)  + 


T  (k)  + 

+ 


R  (k) 


f+(x,k)  =— (k) 


f  ( x ,  -k  )  +  f  (x,k). 


T  (k)  - 


Taking  Wronskians  as  above,  we  obtain: 


(2.6) 


2ik  2ik 


=  If  (x,k),f  (x,k) ] 


T+(k)  T  (k)  '  + 


2ikR  (k) 

— -|k) -  =  (f_(x,k),f+(x,-k)  ] 


2ikR_(k) 

— — =  [f_(x,-k),  f+(x,k)] 


T_^(k)  is  called  the  transmission  coefficient;  R^(k)  are  the  reflection 
coefficients.  The  scattering  matrix,  S(k),  is  given  as: 


(2.7) 


S(k)  = 


T+(k)  R_(k) 


R+(k )  T_(k) 


and  is  well-defined  for  real  k  *  0. 
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The  properties  of  the  scattering  matrix  are  summarized  in  the  following: 


Lemma  2.1.  (c.f.  Theorem  2.1  of  [7))  The  scattering  matrix,  S(k),  has  the 

following  properties: 

(i)  S()c)  is  continuous  for  all  real  k  *  0.  Moreover,  if 

Hull  1  =  f°  ( 1  +  |  x  |  )  2  |u(x  )  |  dx  <  00 
L2 

then  S(k)  is  also  continuous  at  k  =  0. 

(ii)  T+(k)  =  T_(k )  =  T(k ) 

(lii)  TOO  =  T(-k);  R  (k)  =  R^( -k ) 

(iv)  T(k)R  (k)  +  R  (k)T(k)~  =  0; 

| T( k ) |  +  |R+(k)|2  =  | T(k ) | 2  +  | R_ (k ) j  2 

=  1  . 


(v) 


Thus  | T(k ) I ,  | R± (k ) |  <  1. 

Moreover,  S ( k )  is  unitary. 

T(k)  extends  meromorphically  in  the  upper  half-plane  Im  k  >  0. 

T(k)  has  a  finite  number  of  simple  poles  at  k  =  iB  ,  ••♦,i3 

1  n 

where  each  3^  is  a  real,  positive  number;  the  residue  of  T(k) 

at  k  =  i3 .  is 
I 

i[  /  f  (x,iB  )f_(x,iB. )dx]  \ 


-62 , -B2 , • • • ,-32  are  the 
1  2  n 

the  Schrodinger  operator 

away  from  0,iB  ,***,10 
1  n 

at  k=0  as  well). 


eigenvalues  (bound  state  energies)  of 
In  Im  k  >  0,  T(k)  is  continuous 
(if  u  6  ,  then  T(k)  is  continuous 
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( vi  ) 


(vii  ) 


T(k )  =  1  +  0( 1/k) 

as  |k|  +  °°,  Im  k  >  0. 

R±(k)  =  o(1  A) 

as 

|k  j  +  “,  k  real. 

Moreover,  if 

u(x) 

has  £  derivatives  in 

v*> 

k 

•)  as 

1 k |  +  “,  k  real. 

Also,  if  there 

are 

no  eigenvalues. 

2+ 

T(k )  -  1  G 

(Hardy  space) 

and  ! T(k) |  < 

1  for  all  Im  k  >0. 

|  T(k)j  >  0  for 

all 

k  *  0,  Im  k  >  0. 

j  K  |  <  C  |  T(k  )  | 

as 

k  -*■  0  ; 

if  u(x)  e 

then 

ei ther : 

(a)  0  C  C 

1<  1 

T ( k ) | ;  | R  (k ) |  <  C2  <  1 

or 

(b)  T(k) 

=  ak 

+  o(k )  as  k  ♦  0,  Im  1 

R  (k) 

=  -i 

+  Y^k  +  o(k)  as  k  +  i 

Most  of  these  properties  are  direct  consequences  of  (2.5),  (2.6)  and  the 


relation  f^(x,k)  =  f^(x,-k).  The  remainder  come  from  careful  analysis  and 
the  integral  representations  (derived  from  (2.2)  and  (2.5)): 


(  2.8) 


2ikR  (k) 

+  r  -2iky  ,  ,  ,  ,  ,  , 

— 3^-3 —  =  J  e  u(y)m  (y,k)dy 

2ikR  (k) 


-P0 


T(k) 


=  r 


e  2il<;yu(y)m  +  (y,k)dy 


W  =  1  -  2lk  U(y)m±(y'k)dy' 


2ikR  (k) 

Note  that  the  functions  - - -  look  like  Fourier  transforms  of  the 

T(k ) 

potential  u(v),  aside  from  the  factors  m, (y,k).  For  a  full  proof  of  Lemma 


2.  1 ,  see  [ 7 ] . 
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In  [7],  it  is  also  proved  that  m^(x,k)  -  1,  considered  as  functions  of 

k,  belong  to  the  Hardy  space  H  .  The  Fourier  transform  of  m^x.k)  -  1 

with  respect  to  k  plays  an  important  role  in  inverse  scattering  theory. 

As  remarked  in  [8],  the  reflection  coefficient  R+(k)  (or  R_(k))  and 

the  values  8  are  sufficient  to  uniquely  determine  the  entire 

1  n 

scattering  matrix  S(k).  Namely,  since  |T(k)|2  =  1  -  |R+(k)|2  for  k  real, 

given  Rf(k)  we  know  |T(k)|  for  k  real.  Multiplying  by  the  product 
n  k-ifcJ .  n  k-i&. 

II  -i  to  remove  poles,  we  recover  T(k)  *  H  -r  v-g-  by  exponentiating 

j=1  k+i0j  j=1  3 

the  Cauchy  integral  for  its  logarithm.  Then  R_(k)  is  obtained  via  property 
(iv)  in  Lemma  2.1  above. 

Given  R  (k)  and  8.,***,0  ,  does  this  determine  the  potential  u(x)? 

+  in 

It  turns  out  that  n  additional  pieces  of  information  are  needed.  Typically, 

one  specifies  the  so-called  norming  constants,  c j ,  defined  by 

(2.9)  C"1  S  f°  f^(x,i8.  )dx. 

3  _oo  J 

Given  {R  (k ) I k  6  R>  ,  8,, •••,3  ,  and  c,,,,,,c  ,  which  we  call  the 
+  In  in 

scattering  data,  the  potential  u(x)  is  unique  (Levinson's  Theorem).  The 
basic  goal  of  inverse  scattering  theory  is  to  describe  how  one  obtains  the 
potential  from  its  scattering  data. 
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B.  Faddeev's  Inverse  Scattering  Theory 


We  describe  briefly  one  method  of  constructing  the  potential  from  its 

scattering  data.  This  approach,  due  to  Faddeev  [8],  is  based  on  work  of 

Gei'fand  and  Levitan  [10],  Kay  and  Moses  [14],  and  Agranovich  and  Marchenko 

[2).  A  linear  integral  equation  for  the  Fourier  transform  of 

m±(x,k)  -  1,  is  solved  and  u(x)  obtained  by  the  relation: 

8  3 

(2.10)  u(x)  =  -  -5—  3  (x,0+)  =  -5 —  B  (x,  0  ). 

ox  +  ox  — 

More  precisely,  define 

n  -28.  y 

(2.11)  ft  (y)  =  -^  /  R+(k)e21  ydk  +  2  [  c  e  3  . 

j  =  1  3 


The  integral  equation  for  B+(x,y),  y  >  0,  (often  called  the  Gel'fand- 
Levitan-Marchenko  equation)  is: 

(2.12)  B+(x,y)  +  ft+(x+y)  +  /  ft+(x+y+z)B+(x,z)dz  =  0. 

0 

A  similar  equation  for  B_(x,y),  y  <  0,  holds  where 

ft  (y)  H  V  I  R  (k)e-2;Lkydy  +  2  l  c  e  3  ; 


j=1 


(2.13) 


c.  1  =  /°°  f2  (x,i8 .  )dx. 

^  —00  ^ 

We  remark  that  (2.10)  follows  from  the  Fourier  transform  version  of  (2.2), 
while  (2.12)  is  the  Fourier  transform  of  (2.5)  when  suitably  expressed. 

Since  we  will  not  use  the  Faddeev  approach  to  the  inverse  scattering 
problem,  we  do  not  discuss  the  solvability  of  the  Gel ' fand-Levitan-Marchenko 
equation  (2.12)  or  the  properties  of  the  potential  obtained  by  this  method. 

The  interested  reader  is  referred  to  Faddeev  [8]  and  Deift  and  Trubowitz  [7] 
for  more  information  on  this  theory,  and  to  the  works  of  Gardner,  Greene, 
Kruskal,  and  Miura  [9],  Tanaka  [22],  and  Cohen-Murray  [5,6]  for  its 
application  to  the  KdV  initial,  value  problem.  We  merely  note  that  recovery  of 
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the  potential  u(x)  via  (2.10)  does  not  provide  a  simple  means  of  obtaining 
pointwise  estimates  for  u(x)  in  terms  of  R+(k),  as  it  involves  finding 
B+(x,y)  and  differentiating. 

C.  The  Deif t-Trubowitz  Inverse  Scattering  Theory. 

Recently,  P.  Deift  and  E.  Trubowitz  [7]  presented  a  rather  different 
approach  to  the  inverse  scattering  problem.  The  key  to  their  method  is  the 
following  'trace  formula1  (see  also  Newell  [19)) 

li  .00  Oi  If  v  0  p  O 

(2.14)  u(x)  =—  J  k  R  (k)e  m  (x,k)dk  -  4  )  c.3.f  (x,i3). 

tt  +  +  ii+i 

-00  j=1  J  J  J 

Note  that  (2.14)  expresses  the  potential  u(x)  directly  in  terms  of  the 
squares  of  the  Jost  functions,  so  that  the  Schrodinger  equation  may  be  thought 
of  as  a  coupled  system  of  second-order  ordinary  differential  equations  with 
the  singular  boundary  value  m+(x,k)  ~  1  as  x  +  +“>  given.  This  is  the 
basic  idea  in  the  Deif t-Trubowitz  approach.  Carrying  out  this  program 
involves  a  number  of  technical  questions,  which  we  do  not  discuss  in  any 
detail  here.  However,  for  the  reader's  benefit  and  as  an  orientation  for 
Section  5  below,  we  indicate  the  principal  issues. 

First  of  all,  it  is  technically  convenient  to  restrict  attention  to 
potentials  without  bound  states.  In  Section  3  of  [7],  a  spectral  version  of 
Crum's  algorithm  is  used  to  add  or  subtract  bound  states.  Thus  the  problem  of 
recovering  potentials  with  bound  states  from  the  scattering  data  is  in 
principle  reducible  to  the  case  of  potentials  without  bound  states.  The 
method  of  reduction  is  based  on  commutation  ofthe  operator  A*A  where  A  is 

closed.  In  particular,  suppose  Q(x)  is  a  potential  with  bound  state 

..  2  2  2  2  2 
Schrodinger  eigenvalues  -3  <  -3  „<•••<  -3.  •  For  -3  <  -3  let 

n  n-1  1  n 

g(x)  be  a  positive  solution  of 

j  2 

(2.15)  (-  +  Q(X)  +  3  )g(x)  =  0. 

dx 
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Then  it  is  not  hard  to  verify  that,  if  we  define  the  closed  operator  A  as 

(2.16) 


then 


Af  =  <g  h  9"1)f 


*  ,1  2 
A  A  - - -  +  Q(x)  +  3 

dx2 


and 


*  d 

AA  - - -  +  Q(x) 

dx2 

,  2 


„d  ,  ,  ,  02 

2 — -  log  g (x )  +  3 

dx 


+  P(x)  +  3  . 


dx 


In  (7]  it  is  shown  that  A*A  and  AA*  have  the  same  spectrum  except  perhaps 

2 

for  0.  But  -3  ,  which  was  not  an  eigenvalue  for  the  Schrodinger  operator 

with  potential  Q(x),  is  an  eigenvalue  for  the  same  operator  with  potential 
d2  1 

p(x ;  =  Q(x)  -  2  — —  log  g(x)  with  — — -  as  eigenfunction.  This  is  easily 

.2  g  \  x ) 

dx 

seen  since  AA  (— )  -  (-g  — —  g  2  ^ —  g)  —  =  0.  Thus  P(x)  has  an  extra  bound 
g  ^  dx  ^  dx  g 

2 

state  eigenvalue  -3  .  This  describes  how  one  adds  a  bound  state;  reversing 
the  procedure  will  remove  one.  In  this  process,  the  scattering  data  and 
eigenfunctions  are  transformed  in  a  nice  way.  With  Q(’),  P(')  as  above,  the 
transmission  coefficients,  reflection  coefficients,  norming  constants  and 
eigenfunctions  are  related  as  follows: 

'  v*’  ■  CT?  yk> 

V10  ■  -  Mi  V° 

A*  f  (x,i3. ) 


(2.17)  / 


3+3. 


c. 


c. 


j,P  3-3j  j , Q 
c  ,  „  =  2(— J T  ( i 3 ) 

n+1  ,  P  vcr  Q 


where  we  choose  g(x)  above  as 


y  g=t(x,i3)+af(x,i3)  with  a>0. 
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2  2 

Thus  any  -f?  <  -6  and  any  C..,  >  0  are  attainable  this  way.  Moreover,  if 

n  n  +  i 

Q(* )  6  so  is  P(*  )  for  all  m  >  1  and  the  associated  reflection  and 
transmission  coefficients  have  identical  smoothness  and  decay  properties.  In 
[7],  it  is  also  shown  that  n  bound  states  may  be  added  or  removed  in  one 
alqebraic  procedure.  Thus  there  is  no  loss  of  generality  in  considering 
potentials  without  bound  states  in  the  Deif t-Trubowitz  scheme. 

Now  consider  the  Schrodinger  equation 


(2.18) 


m^(x,k)  +  2ikm|(x,k)  =  Q(m+;R)m+(x,k) 
m+(x,k)  ~  1  as  x  ♦  -t* 


where  Q(m+;R)  is  the  right-hand  side  of  (2.14)  (with  no  bound  states). 

.  .  .  _  2ikx  .  , 

n+(x,k)  =  e  m^(x,k). 

We  have  the  equivalent  system: 


Let 


(2.19) 


m^(x,k) 


n^(x,k) 


r.  ~  0 
+ 


-2ikx  ,  ,  . 

e  n+(x,k) 


Q(m+;R)e2i,tXm+(x,k ) 


as  x  -*■  -t“. 


In  order  to  construct  solutions  to  (2.19)  on  the  Banach  space 

B=  j(m,n)|sup  (|m|  +  |n|)  <  “}  a  Lipschitz  estimate  on  the  vector  field  in 
k6R 

(2.19)  is  established.  Then  a  solution  on  semi -inf inite  interval  [M,°°)  is 
constructed,  for  M  large,  by  a  contraction  mapping  argument.  Here  the 
following  condition  on  R+(k)  is  used: 

(2.20)  The  function  F  (y)  =  —  /  R  (k)e2lkydy 

JO  + 

is  absolutely  continuous  with: 

/  |F'  (y)  |  (  1  +  |y|2))dy  <  c(X)  <  00  for  all  <  X  < 

X  + 
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We  remark  that  this  condition  is  rather  natural  when  workinq  with 
potentials,  since  F|(y)  and  u(y)  are  quite  similar  (replace  m+  by  1 
and  ignore  bound  states). 

One  of  the  main  results  of  (7)  is  a  sharp  characterization  theorem  for 
potentials  in  L^.  In  this  class,  necessary  and  sufficient  conditions  for  a 
qiven  set  of  data  to  be  the  scatterinq  data  of  an  potential  are  given  by 

the  following  (Theorem  5.3  of  [7]): 

’sn(k)  s12(k)' 

Theorem  2.2.  A  matrix  S(k)  =  „  „  x  is  a  scatterinq  matrix  a 

1  [S21(k)  S22  j 

real  potential  u(x)  S  L^  without  bound  states  if  and  only  if  conditions  (i)- 
(vii)  of  Lemma  2.1  hold  and  moreover  (2.20)  and  its  analogue  for  R_(k)  (as 
in  (2.13)  above)  hold. 

While  this  theorem  seems  to  suqqest  a  nice  class  of  potentials  for  our 
problem,  the  dispersive  nature  of  the  KdV  equation  leads  to  poor  decay  of 
solutions  as  x  *  -»  for  t  >  0  [1,6).  Therefore,  as  in  Cohen-Murray  [6],  we 
will  reconstruct  the  solution  for  t  >  0  from  x  =  ■¥•>  to  the  left,  usinq  the 
method  sketched  above.  In  this  way,  the  estimate  (2.17)  is  sufficient  and  no 
condition  on  R_(k)  is  needed.  Our  potential  will  be  constructed  for  every 
finite  x  in  this  manner.  One  inherent  difficulty  is  that  the  "left" 
eiqenfunction  f_(x,k,t)  may  not  be  constructed  directly  for  t  >  0  because 
of  the  poor  decay  of  the  potential  as  x  +  -*° .  Thus  the  algebraic  procedure 
for  addinq  or  removinq  bound  states  is  no  longer  available  and  we  are  forced 
to  include  bound  state  contributions  in  the  Deif t-Trubowitz  scheme.  As  we 
shall  see,  this  is  not  a  serious  problem. 


-15- 


3.  The  Inverse  Scattering  Method  for  Solving  the  KdV  Equation:  Evolution 


of  the  Scattering  Data. 

In  a  series  of  papers  culminatinq  in  f9],  Gardner,  Greene,  Kruskal,  and 
Miura  studied  the  KdV  equation 

(*)  qt  +  axxx  -  6qqx  =  0. 

They  discovered  a  remarkable  link  between  this  non-linear  evolution  equation 

and  the  one-parameter  family  of  Schrodinqer  operators 

d2 

(3.1)  L(t)  - - -  +  q(x,t) 

dx 

where  q(x,t),  the  'potential'  in  L(t),  satisfies  the  KdV  equation.  If  we 
consider  the  scattering  data  associated  to  L(t),  namely 

(3.2)  S(t)  =  {R  (k ,  t) :  k  6  R}  u  (S.  (t),c.  (t)}  .  , 

+  3  1  3  =  1 ,*  **  ,n 

then  the  time  evolution  of  the  scattering  data,  assuming  a  priori  that  it  is 

well-defined,  is  linear  and  qiven  by  the  following  well-known  formulae: 

Lemma  3.1.  If  q(x,t)  evolves  according  to  the  KdV  equation  (*),  the 

scattering  data  Sit)  given  in  3.2  satisfy  the  equations: 

f  (i)  ~  R  (k,t)  =  8ik3R  (k,t) 

at  +  + 


(3.3) 


\  (ii) 


tt-  B  .  (t) 
dt  i 


0 


(iii) 


■TT-  C  .  (  t) 
dt  3 


86  3c  .  ( t ) . 
3  3 


Remark.  From  the  discussion  of  the  previous  section  concerning  the 

determination  of  T(k)  and  R  (k)  qiven  R  (k)  and  B.,***,B  ,  it  follows 

-  +  in 


that 


j  a  o 

—  T(k, t)  =  0  and  —  R  (k,t)  =  -8ik  R  (k,t), 
dt  dt  -  - 
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We  present  a  proof  of  Lemma  3.1  due  to  Tanaka  [22],  based  on  the  operator 
formalism  introduced  by  Lax  (16],  which  is  slightly  different  from  the 
oriqinal  derivation  in  [9]. 

As  noted  by  Lax  [10],  the  KdV  equation  is  equivalent  to  the  operator 


equation: 


(3.4) 


—  L(t)  =  [B,L]  =  BL  -  LB 
dt 


where  the  operator  B  is  defined  as  follows  (D  =  -j^-) : 


(3.5) 


B  =  -4D  J  +  3qD  +  3Dq. 


In  the  literature,  the  skew-adjoint  operator  B  is  often  called  the 
Lenard  operator  (after  A.  Lenard)  and  equation  (3.4)  is  often  referred  to  as  a 
Lax  pairinq.  Since  B  is  formally  skew-adjoint,  it  qenerates  a  unitary 
qroup,  so  (3.4)  implies  that  L(t)  and  L(0)  are  unitarily  equivalent  for 
every  t. 

Considering  the  family  of  eigenvalue  problems: 


(3.6) 


L(t)f (x,k,t)  =  k^f (x,k, t) 


and  differentiating  with  respect  to  t,  we  find 


(3.7) 


L(t)[ft  -  Bf]  =  kz[ft-  Bf]. 


If  f(x,k,t)  is  a  Jost  solution  f±(x,k,t),  where  f+(x,k,t)  ~  e  as 
x  +  ¥”  for  every  t  fixed  and  f_(x,k,t)  ~  e  as  x  +  -00  for  every 

t  fixed,  then,  analyzing  the  asymptotic  behavior  of  ffc  -  Bf,  we  find: 


(  3.8) 


(f+)t  -  Bf+  =  -4ik3f+ 


(f_)t  -  Bf _  =  4ikJf _ 


Here  we  use  the  assumption  q(x,t),  q  (x,t)  -*■  0  as  |x|  +  °°,  t  fixed. 
Differentiating  the  relation  (2.5)  with  respect  to  t  and  usinq  (3.8),  we 
obtain  (3.3) (i).  (Alternatively,  one  can  differentiate  the  Wronskians  in 


(2.6)  and  obtain  this  result.)  By  our  previous  remark,  L(t)  is  unitarily 


-17- 


equivalent  to  L(0),  thus  3  ^  ( t )  =  3  1 0 ) .  To  obtain  (3.3)(iii),  we 

differentiate  f  (x,i8.,t)  =a.(t)f  (x,i8.,t)  and  use  (3.8). 

+  3  3-3 

An  interpretation  of  the  linear  flow  (3.3)  in  the  context  of  completely 
inteqrable  Hamiltonian  systems  was  given  by  Zakharov  and  Faddeev  [28]. 
Introducing  a  Hamiltonian  structure,  they  show  that  the  variables 
[  P(k )  =  -  ^  «n(1  -  jR_(k)|2) 


R  (k) 

900  = 


(3.9) 


P.  =  6 
3  3 


qj  2en (lCj  dk  ^T(k)^  ^k=i8  ) 


form  action-anqle  variables  for  the  KdV  Hamiltonian.  In  fact,  they  derive  the 
expression 


(3.10) 


H[u]  =8  r  k3P(k)dk  -  -p -  l  p?5/2 


which  is  equivalent  to  (3.3)  reexpressed  in  terms  of  F  (k),  c.,  8.  (the 

3  3 

alternative  set  of  scattering  data  using  normalizations  at  x  =  -<x>). 
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4.  Recovering  the  Potential  at  Later  Times. 

Shifting  our  viewpoint  somewhat,  we  now  consider  the  linear  flow  on  the 
scattering  data  induced  by  the  nonlinear  KdV  equation.  We  show  that  under 
suitable  restrictions,  the  evolved  scattering  data  is  sufficiently  "well- 
behaved"  to  permit  recovery  of  an  associated  potential  by  the  inverse 
scattering  methods  of  Section  2  above.  To  this  end,  we  make  the  following 
def ini tion: 


(4.1)  A  function  q(x,t)  is  called  a  generalized  solution  of  the  KdV  equation 

(in  the  sense  of  inverse  scattering  theory)  for  all  t  >  0  if  it  is  the 

potential  corresponding  to  the  scattering  data: 

..3  #.t 

(4.2)  S (t)  =  {R+(k,0)e  |k  6  r}  u  {B^,c^(0)e  3  |j  =  1,***,n} 

for  every  t  >  0. 

To  avoid  confusion,  the  term  "weak  solution"  will  be  used  for  solutions 
in  the  sense  of  distributions  and  "generalized  solution"  will  always  refer  to 
the  notion  (4.1).  The  chief  result  of  this  section  is  the  following  theorem. 
Theorem  4.1.  Given  Q(x)  G  L^,  there  exists  a  unique  generalized  solution 
(in  the  sense  of  inverse  scattering)  for  t  >  0  to  the  KdV  equation 


(*) 


qt  +  W  -  =  0 


with  scattering  data  S( t)  as  in  (4.2)  where  the  initial  values,  S(0),  are 
the  scattering  data  for  Q(x). 

The  proof  of  the  theorem  consists  of  three  steps,  the  first  of  which  (our 
Lemma  4.2)  appears  in  (7]  (Theorem  4.7  and  Remark  4.5  thereafter). 

Lemma  4.2.  If  Q(x)  e  L^,  then  the  function  H(x)  defined  by 


(4.  3) 


H(x)  i  CgM  /3  R+(k,D)e2lkXdk 

-a 


i  s 


also  in  L  ,  where  RL(k,0)  is  the  reflection  coefficient  for  p(x) 
a  + 


and  C^M  denotes  the  Cesaro  mean. 
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r  v/V  " 


For  the  sake  of  completeness,  we  sketch  the  proof  appearing  in  [7], 
Use  the  following  version  of  the  trace  formula: 


(4.4) 


Q(x )  =  H(x)  +2  /*  H ( x  vy ) B ( x , y ) dy  +  /"  H(x+y) (B*B) (x,y)dy 
0  0 

n  -26.x 

-  1  4  c .  6.e  3  in  (x,i6 .  ) 

j=1  3  3  +  3 

where  B(x,y)  (resp.  B*B)  is  the  inverse  Fourier  transform  (in  k)  of 
m+(x,k)  -  1  (resp.  (m+(x,k)-1 ) 2) ,  and  6j ,c  are  the  jth  bound  state  and 
norming  constant  for  Q(x).  Recall  that,  as  in  Section  2,  Q(x)  determines 
m+(x,k ) . 

Solving  (4.4)  by  iteration  for  H(x)  leads  directly  to  an  estimate  of 
the  form 

(4.5)  «H(*)«  <  C(  a)  #Q(  • )  II 

L1  L 

a  a 


with  bounded  constants  C(a)  for  all  a  >  1.  Thus  Q( • )  6  implies 

HO  e  l|. 

Remark .  The  trace  formula  (4.4)  and  its  time-dependent  generalization  will 

play  an  important  role  in  the  discussion  of  Section  5,  in  which  further 

properties  of  the  solution  constructed  in  Theorem  4.1  are  deduced. 

The  second  step  in  the  proof  of  Theorem  4.1  concerns  the  behavior  of  the 

Sik  ^  t 

Fourier  transform  of  2ikR+(k,0)  •  e  for  t  >  0  and  is  crucial  to  the 

analysis  of  Section  5  below. 

Lemma  4.3.  The  function  h(x,t)  defined  by 


(4.5) 


h(x,t)  =  —  2ikR  (k,0)e 


2ikx+8ik  t 


dk 


is  well-defined  for  t  >  0  and  satisfies: 
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h  ( *  ,  t )  6  for  every  t  >  0  fixed  and  every  6  >  0 


with  t  +  h(*,t)  continuous  in  this  topology 


t  +  h  (•  ,  t)  6  C(]0,«[,X.) 

4 


In  fact,  lim  (1  +  Ixjj^fxjt)  =  0  for  t  >  0  fixed 


(iii)  h(‘ ,t)  6  L 


r  s 

All  of  the  above  hold  for  3  fc3  h(x,t)  provided  we 
subtract  V2  ( 3r  +s )  from  every  lower  index  (the  poly¬ 
nomial  weighting)  i.e. 


3*3Sh(x,t)  e  x1 
tX  1 3/4-5  -2|±i 


for  every  t  >  0  fixed 


h(xft)  +  H(x)  in  X2-&  aS  t  *  0 


h(x,t)=  0(t  )  as  t  +  4°°  for  fixed  x. 


(vii  )  h(ct  +  £  ,  t)  = 


o(t"13/3)  if  c  >  0 
0(t"  ^  )  if  c  <  0 


as  t  +  -f30  (£  fixed) . 


Remark.  (iv)  implies  that  for  3r+s  <  5, 


3t3x  is  continuous. 

Thus  h(x,t)  is,  for  t  >  0,  a  classical  solution  of  the  linear  equation: 


(4.7) 


ht  +  hxxx  =  0 


and  (v)  describes  the  sense  in  which  the  initial  data  is  taken  on. 


-21- 


the 


Proof  of  Lemma  4.3.  We  first  note  that  for  potentials  in  L^  , 

-  4 

associated  reflection  coefficient  R+(k,0)  is  C2  (see  [7]).  Given 
and  x,  both  fixed  and  finite,  we  may  pick  a  finite  kg  with 

kQ  >  ( ^  »  which  implies  x  +  12k2t  >  0  for  |k|  >  k^. 

Inteqrate  (4.5)  by  parts  in  |k|  >  k  ,  brinqing  down  a  factor  of  x  + 
in  the  denominator,  to  show  that  the  inteqral  in  (4.5)  is  well-defined 
any  t  >  0  and  x  fixed. 

In  order  to  verify  the  properties  of  h(x,t),  we  re-express  (4.5), 
(4.3)  above  to  replace  2ikR+(k,0)  by  the  inverse  Fourier  transform  of 
H( x ) .  We  then  have: 

,  .  ,  ..  1  f“  %  2ik(x-y  )+8ik3t 

(4.8)  h(x, t)  =—  J  J  H(y)e  7  dydk 

_po  _po 

where  the  y-inteural  is  absolutely  convergent  and,  as  remarked  above, 
iterated  integral  converges.  Thus  h(x,t)  is  the  convolution  of  H(x) 
the  fundamental  solution 


..  1  f”  2ik(x-y)+8ik3t ,, 

(4.9)  E(x-y,t)  =  —  J  e  7  dk 


2iK(^W)  .  I1"3 


-  r 


(  3t ) 


1/3'  3 


dK 


-  ( 3t )”* '3  ) 


where  Ai(z)  is  the  well-known  Airy  function. 
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t  >  0 


1  2k2t 
for 

using 


the 

with 


r 


We  shall  make  use  of  the  following  properties  of  Ai(z)  (c.f.  [20]  for 


example ) : 


( a )  Ai  "  ( z  )  =  z  Ai  ( z ) 

(b)  | Ai ( z  )  |  <  1  for  all  real  z 

(c)  For  z  >  0,  Ai(z)  is  a  non-negative,  monotone  decreasing 
function  of  z  with  asymptotic  rate  of  decay  proportional  to 

-  V4  .  while  Ai '  decays  like  z+  ^4  -5  where 


C  S  2/3  z3/2 

(d)  As  z  ♦  ,  Ai(z)  is  hiqhly  oscillatory  with  an  algebraically 

decaying  envelope.  In  fact, 

Ai(z)  ~  — Vr  |z|  ^  {  cos  (?  -  it/4)  +  0(— )} 
tt  /2  l 


where  C  =  2/3  | z 


3/2 


and  Ai ' (z  ) 


V4 


1  {  sin  (C  -  if/4)  +  0(^-)J. 

C 


„  V2 


Using  (a),  corresponding  estimates  on  higher  derivatives  of  Ai  hold. 
We  now  show  that  (4.6)  holds: 


r  a  ♦  i*i) 


1 3/ 4  — <5 


h(x, t) |dx 


(4.11) 


=  /  (1  +  |x|  )  1  3/4-6  |1  1  Ai(iS_y)H(y)dy|dx 


1  /3 

where  €  =  ( 3 t )  >  0,  fixed 


f n  +  |x| 

,i 3/4-6 1 

1  e  /  Al(  £  ) H(y )dy |  +  |e  /  Ai(  £  ,H(y)dv|}dx 

a 

r<  1  +  ixi 

j  1  3/4-6  { 

'  -CO  X 

|  h 1 ( x , t ) |  +  |h2 (x, t ) | } dx. 

a 
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Now  we  have 


(4.12)  e  •  /"*  (1  ♦  |x|  )°|h1  (x,t)|dx 
a 

<  /  (1  +  |x|)°  /a  |  Ai  (-~~)n(y )  |  dydx  +  /  (1  +  |  x  | )  a/X  |  Ai  (-^L}H(y )  |dydx 
a  a  a 

=  /ady  f°  dx  |Ai(— |  |  H(y)|(1  +  |x|)“  +  /“ dy  /“dx  |Ai(^-)  |  |  H(y )  |  ( 1  +  |x  |  )  ° 
-“a  ay' 


<  /a | H( y ) | dy  •  /“( 1  +  | x  j ) °Ai }dx  +  /  ( 1+|y | ) a| H(y ) |dt  /  ( 1+|z | ) °Ai (z/e)dz 
-“a  a  0 


< 


We  remark  that  these  estimates  give  a  bound  on  llh  for  all 

a(a ) 

a  <  4  uniformly  in  0  <  t  <  T  <  »  by  rescaling. 

Similarly, 

(4.13)  /*(  1  +  |x  |  )“|h2(x,t)  |dx  =  /"( 1+|x  |  )“  /“|Ai(^-)  |  |  H(y )  |dydx 

a  ax 


=  /“  /°( 1  + 1 x  |  )a|Ai(z)  |  |H(x-ez)  |dzdx 
a 


<  r /°<hx  -ez | )a | Ai (z ) | | H(x-ez ) | dzdx  (assuming  w.A.o.g.  a  >  0) 


<  /°} Ai ( z )  |  /"  <1+|s|)°|H(C)|dCdz 

-00  a-ez 


4  J°  |  Ai  ( z )  |  •  (1  +  |a-ez|)a  4dz  •  »H(,)II 

L^(a) 

4 


<  00  by  the  decay  rate  on  Ai  if  a  <  13/4. 
Note  that  this  is  not  uniform  as  e  +  0. 
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To  estimate  (1  +  we  proceed  as  follows  (0  <  6  <  1) : 


(4.14)  ( 1+|x| )4h(x,t)  =  (1+|x|)4  /x^ 1 Ai(^^-)H(y  )dy 


.  I1.|x|  -  1  »i(ifr)H(y)dy. 


x( 1-6 ) 


The  second  term  goes  to  0  as  x  >  +“,  since  (1+|y|)4H(y)  G  L1  and 


Ai  6  L  .  In  the  first  term,  y  <  (1-6)x  so 


x-y  >  6x 


By  monotonicity, 


/  X -V  ^  (  6  X 

Ai(— —  J  <  Ai(_— J,  which  decays  exponentially  as  x  ♦  +00.  Thus 
lim  ( 1 +|x | )4h(x,t)  =  0  for  every  t  >  0. 

x++.<» 

To  obtain  decay  of  h(x,t)  as  x  +  -“,  we  make  a  similar  estimate: 


(4.15)  (  1  +  |  x  |  )  ^4|h(x,t)|  =  (1+|x|)  ^  |  /X(1  Ai  H(  y  )dy  | 

—00 

+  { 1  + 1 x  I  )  1/4  I  r  -  Ai (— jp^)H(y  )dy. 
x{ 1-6) 

As  x  +  the  first  term  is  0( < 1  +  lx  I )-1 3/4)  since  H(*)  G  L1.  while  for 

4 

the  latter  term  (1+|xj)  ^  |  Ai(X£^-)  j  is  bounded  and  H(')  G  L,1.  Hence 
( 1 + | x | )  |h (x, t ) |  G  L  for  every  t  >  0. 

The  extension  of  these  estimates  to  x  and  t  derivatives,  as  well  as 
their  continuity  in  t  for  t  >  0,  follows  from  the  analyticity  and 
estimates  on  the  Airy  function. 


Asymptotics  in  t  are  derived  as  follows: 


(4.16) 


(  3t)  1//3h(x,t)  =  j  Ai  ( — X  ^,rr)H(y)dy. 

(3t)VJ 
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Note  that  this  aqrees  with  the  alternate  method  of  estimate,  namely  stationary 
phase  analysis  of  (4.5). 
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For  c  >  0,  we  have: 


(4.18)  I  (3t)1/3h(ct-^  )  I  <  /  2  | Ai |  |l-f(y)|dy 

-*>  ( 3t) 


+  !“  jAi(~— ~yr-)  j  |h(v)  |dy 

ct  r  ( 3t )  ' 


<  sup  |  Ai( - ^TTr)  I  /  |  H(y )  [ dy 

.ct  -*> 

Z>— 


r  n  +  jy|)4|H(y)|dy  =  o(t  ), 


, .  i ct  r | .4  ct  r 
d+!_— £|)  _4C 

-13/3 

hence  we  have  in  this  case  h(ct+£,t)  =  oft  )  as  t  + 

To  show  that  h(x,t)  +  H(x)  in  as  t  +  0,  first  we  consider  (as 

usual  £  =  (3t)1/3) 


(4.19) 


/e  Ai (y) [H(x-£y)  -  H(x)]dy. 


Now  H (• )  6  implies  that,  since  translation  is  continuous. 


H(x-ey)  -  H ( x ) II  1  =  o ( H e yll  )  as  jey[  +  0. 


On  the  criven  interval,  | E v  I  <  e  2  while  the  interval  lenqth  is  2£ 
Thus  the  expression  (4.19)  is  o(1)  as  e  *  0  as  an  element  in  L^. 
Now  we  have: 


(4.20) 


-r1 /2 

h(x,t)  =  J  Ai (y )  H( x-e y)dy 


+  /E  Ai (y)H(x-ey )dy  +  /  Ai (y )H(x-£y )dy. 


-2  7- 


n 


The  first  term  is  estimated  in  ^ ( a )  by: 


_  -L 1/2 

(4.21)  f°  /  E  |Ai(y) | (1+|x| )  3~  | H(x-ey ) | dydx 

a  -® 

"  -V2 

<  c|e) 1/<8  •  /“  /  e  (l  +  jx|)3-  |H(x-£y)  | dydx 


<  c|e|1//8  /** (1  +  |xj)3"6  /  1  |H(z)|dzdx 
a  x+e  '2 


C I  C  |  1/8  /”  (1+|xj)  1  ^dx  •  IH(  •  )  N  1 


I*' (a) 
4 


=  0(|£|V«). 


/*  /"  |Ai  (y )  I  ( 1  +  |x|  )3_l5 1  H(x-ey)  jdydx 

a  rV 2 

<  /"  /"  |Ai(y)|(1+|x-ey|)3_6(l+|ey|)3_6|H(x-ey)|dydx 

"  7 1/2 

<  /”  | Ai (y ) | (1 +| ey | ) 3  6dy 

7  V2 


which  tends  to  0  exponentially  fast  as  £  +  0. 


h(x,t)  +  lim  /£  Ai(y)H(x)dy  =  H(x) 

-  iv. 


in  X'  .  as  t  +  0. 

3-o 

This  completes  the  proof  of  Lemma  4.3. 
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Remark.  If  we  assume  more  generally  Q(x)  C  L^,  then  H(x)  t>  and  the 
corresponding  results  for  h(x,t)  are: 

(4.21)  (i)  hC.t)  x-_3/4_6  n  x”  n  l”/4 

for  t  >  0  fixed,  with  t  h(*,t)  continuous  for  t  >  0. 


for 

c  <  0 

(ii) 

for 

c  >  0 

(iii)  h(*,t)  *  H(  • )  in  X*  x  as  t  +  +“>  (£  fixed). 

a- l-o 

i  °» 

In  particular,  if  Q  6  L  with  compact  support,  then  h(x,t)  is  C  for 

t  >  0,  h  decays  faster  than  any  power  of  x  as  x  +*•  for  fixed  t  and 

h  decays  faster  than  any  power  of  t  in  moving  frames  with  c  >  0  as 

t  +  +"  (c.f.  [6]).  Better  decay  as  x  +  -°°  and  in  moving  frames  for 

c  <  0  is  related  to  additional  smoothness  of  Q(x)  and  hence  H(x). 

Using  Lemma  4.3  and  the  results  of  [7]  (c.f.  (2.17)  above),  it  follows 

that  for  every  t  >  0,  a  potential  q(x,t)  may  be  recovered  via  inverse 

scattering  starting  from  +00  and  solving  to  the  left  for  every  finite  x. 

Thus  a  generalized  solution  (in  the  sense  of  inverse  scattering)  exists  and, 

for  each  fixed  t  >  0,  q(*,t)  6  Xj  ^  In  ttle  next  section,  we  strengthen 

this  result  to  include  smoothness  and  decay  of  the  generalized  solution. 


5.  Regularity  and  Decay  of  the  Generalized  Solution. 

By  analyzing  the  inverse  scattering  method  of  Deift  and  Trubowitz  [7] 
more  closely,  smoothness  and  spatial  decay  properties  of  the  generalized 
solution  constructed  above  will  he  obtained.  The  solution  behaves  quite 
similarly  to  h(x,t),  as  discussed  in  Lemma  4.3  above,  except  for  the  presence 
of  a  finite  sum  of  "soli ton”  terms.  We  shall  prove  the  following  theorem. 
Theorem  5.1.  The  generalized  solution  q(x,t)  constructed  above  satisfies: 
(5.1)  fr(i)  q  (*  / 1)  S  xJ3/4-6  n  x“  for  t  >  0 

and  t  ♦  q(*,t)  is  continuous  for  t  >  0 


(ii) 


3  *3  8  q  ( •  ,  t )  e  x1 

t  X  , 


,  n  X 

r  (S+3ri  .  r  3r+s-i 

3/4-6  _( — _ — J  4-1— 2~J 

(with  continuity  in  t) 

1 

(iii)  q ( x , t )  -  0(x)  in  X^. 


Our  proof  is  essentially  a  parameterized  (by  t)  version  of  the  Deift- 
Trubowitz  approach  to  inverse  scattering  [7],  sketched  in  Section  2  above.  We 
begin  by  rewriting  the  time-dependent  analogues  of  the  trace  formula  (2.14) 
and  the  Volterra  equation  for  m+  (2.2)  in  terms  of  b(x,y,t),  which  is  the 
inverse  Fourier  transform  of  ir  (x,k,t)-1  with  respect  to  k.  Since  m+-1 
helonqs  to  the  Hardy  space  H2  ,  b(x,y,t)  =  0  for  y  <  0.  We  have  the 
Volterra  eguation 

(5.2)  h ( x, y , t )  =  I"  q(z,t)dz  +  /ydz  /"  dw  q(w,t)b(w,z,t),  y>  0 

x+v  0  x+y-z 


and  the  "trace  formula 


(5.3)  q(x,t)  =  h(x,t)  +  2  /  h(x+y,t)b(x,y,t)dy 

0 


+  /°°  h(x+y,t)(b*b)(x,y,t)dy 
0 

n  -20.x+8e3t 

h(x,t)  =  h(x,t)  -  'i  40.c.(O)e  3 

j=1  3  3 


where,  as  in  Section  4  above,  h(x,t)  is  the  Fourier  transform  of 
2ikR+(k, 0)e8ik  (see  (4.5)  above).  Note  that  the  addition  of  the 
exponential  terms  does  not  alter  the  results  of  Lemmas  4.2,  4.3.  We  regard 
(5.2),  (5.3)  as  the  fundamental  equations  to  be  solved,  and  as  in  (7),  we 
solve  them  by  iteration  to  the  right  of  a  sufficiently  large  point  (depending 
on  h(x,t))  and  then  argue  that  the  solution  continues  to  exist  for  all 
finite  x  values,  by  using  a  local  Lipschitz  estimate. 

One  important  step  in  this  procedure  is  to  obtain  estimates  on  the 
solution  of  (5.2)  if  we  assume  q(x,t)  is  known.  We  introduce  auxiliary 
functions  n(x,t),  Y(x,t)  assuming  q(*,t)  6  X3  for  every  t  >  0,  as 
follows : 


(5.4) 


n(x,t)  =  /  |q(y,t) |dy 


Y (x,t)  =  f°(y-x) |q(y,t) |dy. 


Note  that  n,Y  are  non-increasing  functions  of  x  for  each  fixed  t  which 


are  finite  for  all  x  >  -°°  and  that 

/  n(y,t)dy  =  /  dy  /°°  dz  |q(z,t)| 
x  x  y 

(5.5) 


=  /  dz |q(z, t ) | (z-x )  =  Y(x,t). 
x 


Similarly,  if 


3lq(' 


,t)  e  x 


i ' 


we  may  introduce 
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(5.6)  f  X(x,t)  i  J  |yjr(y,t)|dy 


L  u(x,t)  i  f* (y-x) ||^(y,t) |dy. 
x 

With  these  functions  thus  defined,  we  prove  the  analogue  of  Lemma  2.3  of  [7). 
Lemma  5.2.  If  q(*,t)  6  X^  for  every  t  >  0  and  ■^■(•ft)  6  *•]  for  every 
t  >  0,  the  integral  equation  (5.2)  has  a  unique  solution  b(x,y,t) 


satisfying: 
(5.7)  f  (i) 


|b(x,y,t)|  <  e^X,t*n(x+y,t) 


(ii)  b(x,y,t)  is  absolutely  continuous  in  x  and  y  with 

||p-(x,y,t)  +  q(x+y,t)|  <  eY(x't)n(x,t)  •  n(x+y,t) 

||^-(x,y,t)  +  q(x+y,t)|  <  2eY^x,t^n(x,t)n(x  +  y,t) 

(iii)  b  solves  the  'wave  equation' 

J  3^  g2 

\  — ^(x.y.t)  -  b(x,y,t)  =  q(x,t)b(x,y,t) 

3x^ 

with  -  b(x,0+, t)  =  -  b(x,0+,t)  =q(x,t) 

(iv)  ||^(x,y,t)|  <  e^X,t)U(x+y,t)  +  eYtx,t)  •  n(x+y,t)y(x,t)} 

(v)  If  (•|^-J^q(x, t)  6  X^  for  every  t  >  0  and  0  <  j  <  £ 

then  b  has  mixed  deriatives  in  x  and  y  up  to  order 
1+1  obeying  estimates  analogous  to  (ii)  above. 

Proof .  (i),  (ii)  and  (iii)  are  proved  in  Lemma  2.3  of  [7],  so  we  merely 

sketch  their  proof. 


Solving  (5.2)  by  iteration  with 
00 

(5.8)  [  b(x,y,t)  =  l  b . (x,y,t) 

3=0  3 


b  (x,y,t)  =  /  q(z, t)dz 
x+y 


j.  (x,y,t)  =  fy  dz  /  dw  q(w,t)b. (w,z,t) 
•5+  0  x+y-z  3 


leads  to  the  inductive  estimate  [7} 


I b j  ( x ,  y ,  t )  |  <  ~'Y|' j  ~  n(x+y,t) 


which  proves  (i). 


Differentiating  (5.2)  and  using  the  above  estimate,  we  have 


(5.10) 


r 


1 3^-  b(x,y,t)  +  q(x+y,t) 


=  |  /y  dz  q (x+y-z , t)b (x+y-z, z, t) j 
0 

*  /y  dz |q(x+y-z, t)  |eY^x+y-z't'n(x+y,t) 
0 

4  eY<X,t,n(x+y,t)  /y|q(x+y-z,t) | 

0 

<  eY(x't)n(x+y,t)n(x,t) 


with  a  similar  estimate  for  the  y  derivative. 

(iii)  is  a  direct  calculation. 

To  establish  (iv),  differentiate  (5.2)  with  respect  to  t  ,  to  obtain 
3b  ,  .  r°°  3q 


(5.11) 


3t 


(x,y,t)  =  f  ^  (z,t)dz 


x+y 


+  r dz  /  dw  (|^-(w,z,t)q(w,t)  +  b(w,z,t)|5-(w,t)  ). 


„  at 

0  x+y-z 

Again  we  solve  by  iteration,  writing 

r 

(5.12) 


at 


b  (x , y , t )  =  l  £ . (x,y,t) 
j=0  J 


r *»  3q 

C  (x,y,t)  =  J  dz  (z,t) 
x  +y 


+  /Y  dz  /  dw  b(w,z,t)  (w,t) 


0  x+y-z 


C 


3+1 


(x ,y , t ) 


V 


/ydz  /  dw  q(w, t)C . (w,z, t). 
0  x+y-z  3 
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Then  |c_(x,y,t)|  4  A(x+y,t)  +  J"ydz  J*  dw|-^.(w,t)  n(w+z,t) 

0  x+y-z 

<  A(x+y,t)  +  n(x+y,t)u(x,t)  •  e^x,t/ 

and 

lCn(X'y't>!  <  l5oU'y't)l  *  ••Y('n1t)) 

(established  just  as  before),  which  proves  (iv).  (v)  is  proved  inductively 

by  applying  the  same  ideas  as  in  (ii)  to  higher  derivatives.  Namely,  from  the 
relation 

(5.13)  (|^)j+1b(x,y,t)  +  (|^Pq<x+y,t) 

=  -  /y(|^-)-’ (q(x+y-z,  t)b(x+y-z,z,  t)  )dz 
0  x 

and  estimates  on  derivatives  of  q  and  b  of  lower  order,  we  can  readily 

estimate  (-^)^+1b(x,y,t).  This  completes  the  proof  of  Lemma  5.2. 

We  shall  also  require  estimates  on  the  difference  between  the  functions 

b^)(x,y,t),  j  =  1,2,  corresponding  to  two  potentials  q^,q*2^.  Define 

5n(x,t)  i  /**|q(1  )  (y,t)  -  qt2)(y*Z)|dy 
x 

(5.14)  6y(x,t)  =/  (y-x) |q( 1  * (y,t)  -  q(2> (y,t) jdy 

X 

^  a  (1)  a„(2) 

<SA(x,t)  =  r  |-  ■q-3-t--(y,t) - g^—  (y,t)|dy 

X 

^  3d)  3(2) 

6m (x , t )  i  J  (y-x )  |— — (y,t) - ^ — (Y»t)|dy 

x 

for  notational  ease.  We  have: 

Lemma  5.3.  (c.f.  Lemma  2.4  of  (7))  If  q^^(x,t),  q^2^(x,t)  are  two 

potentials  as  in  Lemma  5.2  above,  then  the  corresponding  solutions 
b^Nx,y,t),  b^2*(x,y,t)  of  (5.2)  satisfy: 
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have 


(5.16) 


P(x,y,t)  =  b(1)(x,y,t)  -  b(2)(x,y,t) 

r®°  ( 1  )  ( 2 ) 

=  I  (q'  1Z/t)  -  q'  '(zft))dz 

x+y 


+  /ydz  /**  dw{q  ^ 1  *  (w,  t)  (b  ^ 1  ^  (w,  z,  t)  -  b^2^(w,z,t)) 

0  x+y-z 

+  (q  ^ 1 ^ ( w, t)  -  q^2' (w, t) )b<2' (w,z,t)} . 

00 

Solvinq  for  p  =  £  p.(x,y,t)  leads  to  the  usual  inductive  estimate 

j=0  D 

(Y  .  ( x,  t) ) ^ 

(5.17)  |p_.(x,y,t)|  <  |p  Q(x,  y,  t )  |  •  - — - 

where 

(5.18)  P0(x,y,t)  =  f°°  (o  ^ 1  ^  ( z ,  t )  -  qt2'(z,t))dz 

x+y 

+  /ydz  /  dw(q(1)(w,t)  -  q ( 2) (w, t) )b( 2> (w,z, t). 

0  x+y-z 

Thus 

|p  Q(x,y, t) |  <  f  |q(1)(z,t)  -  q^2)(z,t)|dz 
x+y 


+  /ydz  f°  dw |g^  ^  ( w,  t ) 
0  x+y-z 


Y  (v,t) 

q '  (w,t) |e  n2<w+z,t). 


(i)  follows  by  interchanging  the  order  in  the  double  integral  and  integrating 
n^(w+z,t)  while  (ii)  comes  from  using  the  monotonicity  of  n  which  yields 
the  n2(x+y,t)  factor. 

(iii)  and  (iv)  are  derived  similarly.  Differentiating  (5.16)  with 
respect  to  t  gives 


-36- 


(5.19) 


3b 


(  1  ) 


(x'y/t)  — (x,y,t)  =  (-fi — (z,t)  -  — (z,t)j 


(2) 


Tt 


Tt 


x+y 


3q 

3t 


(1) 


9q 

at 


(2) 


+  /ydz  f  dwJ 
0  x+y-z  ' 


^  ( 1 )  to) 

- jj"£ — (w,t)[b  (w,z,t)  -  bV  ( w,  z  ,  t ) ) 


(1)  3b(1>  3b(2) 

+  q  (w,t)  [— — (w,z,t) - ^ — (w, z, t ) ] 


+  — (wft) 


(D 


3  (2) 

(w,t))bl2)(w,z,t) 


(2) 


+  (q  (w,t)  -  ql  (w,t))  — — (w,z,t)}. 


3b(1)  3b(2) 


Solving  iteratively  for  — —  again  leads  to  an  inductive  estimate 
involving  (Y  ^  (x ,  t ) ) -1 /j  1  ,  but  now  the  leading  term  is 


r, 3q(1)  3q(2) 

/  l"y£ — (Z,t) - yj; — (z,  t)  )dz 

x+y 


+  Jy dz  /  dw / 
0  x+y-z 


-yy — (w, t) [b^ 1  * (w, z, t)  -  b(2* (wrz,t) ] 


+  ^~Tt - (w,t*  ~  —$t - IW'  2)  (w,z,  t) 


v. 


(  2  1 

+  (q(1>(w,t)  -  q(2)(wft))-^ — (w,z,t) 


Using  (l),  (ii)  and  the  estimates  of  Lemma  5.2,  depending  on  whether  we 
integrate  out  the  bounds  for  the  b's  or  take  their  suprema  over  z  ,  we 
obtain  (in )  and  (iv)  above.  This  proves  the  lemma. 
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The  above  lemmas  will  be  used  to  show,  in  a  manner  analogous  to  17],  that 
the  mapping 

(5.20)  q  |+  *(q) 

given  by  finding  b  from  q  via  (5.2)  and  then  defining  $(q)  from  b  via 

(5.3)  is  a  contraction  on  a  suitable  ball  in  a  Banach  space.  In  [7],  this  space 

is  {Q(x)  :  /  | x | | Q( x ) | dx  <  l)  where  a  is  chosen  so  that  /  |x||H(x)|dx  is 
a  a 

sufficiently  small  (recall  H(x)  is  the  inverse  Fourier  transform  of 

2ikR+(k,0)).  Since  we  are  interested  in  smoothness  and  decay  properties  of 

our  generalized  solution,  we  shall  include  these  in  our  Banach  space. 

Before  discussing  the  mapping  q  )■►  ®(q)  more  fully,  we  require  some 

estimates  on  the  relation  (5.3),  where  we  regard  b(x,y,t)  as  known. 

Lemma  5.4.  Suppose  functions  Y(x,t),  h(x,t),  X(x,t),  u(x,t)  exist 

satisfying  (5.5)  and  its  analogue  for  X,u  and  such  that  the  estimates  (5.7)  (i) 

and  (iv)  hold.  Then  if  q(x,t)  is  given  by  (5.3),  we  have,  for  every  •<  >  0 


lq(*  »t) «  <  (1  +  Y(M,t)eY(M,t))2  •  llh(  *,t)  I 

La(M)  La(M) 


Uq  ( - ,  t )  U  „  <  (1  +  Y(M,t)eY(M,t)}2  Ilh  (  • ,  t )  II 


La(M) 


La(M) 


■|f(*rt)l  1  4  (1+Y(M,t)eYlM,t^  ) ^  {»£(., t)« 

3t  L>)  3t  L>) 


.  Y (M, t ) .2  ru3h, 


+  2U(M,t)eY(M,t)«hC,t)U  1  j 


La(M) 


'Ifc.t)!  1  <  (l  +  Y(M,t)eY(M,t)]2  l||c,t)l  1 


L  (M ) 
a 


La(M) 


+  2Y(M,t)( 1+Y(M,t)eYCM,t) ) ( 1+n(M,t)eY<M,t) Uh( *,t) a 


L  (M) 
a 
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Proof  of  Lemma  5.4.  From  (5.3), 

q(x,t)  =  h(x,t)  +  2  /  h(x+y,t)b(x,y,t)dy  +  /  h(x+y,t) (b*b) (x, 

0  0 

Thus 

(5.22)  /”( 1+|x| )a|q(x,t) |dx  <  /  ( 1 +| x | ) a|h (x, t ) |dx 

M  M 

+  /“( 1+jx| )a|2  /“  h(x+y,t)b(x,y,t)dy |dx 
M  0 


+  /“(1+|x|a)|  /  h (x+y, t) (b*b ) (x,y, t )dy |dx 
M  0 


<  Ilh  ( • ,  t)  I  +  2  /“(1+jx|  )“  /"jh(x+y,t)  |eY(x,t)n(x+y,t)dydx 

L  (M)  M  0 

a 


+  /"(1+|x|)a  /*|h(x+y,t)| |b*b(x,y,t)|dydx 
M  0 


<  lh(*,t )»  +  2  /"n(x,t)eY(x,t)  /“(7  +  |x+y|)a|h(x+y,t)|dydx 

I/(M)  M  0 

+  r Y(x,t)n(x,t)e2Y(X,t>  /“>(1  +  |x+y|)°‘|h(x+y,t)|dydx 


«;  Ilh  (  •  ,  t )  1  ,  ll  +  2Y  (M,  t  )eY  +  Y2(M,t)e2Y(M,t)}. 


La(M) 


Similarly, 


(5.23  II  q  (  *  ,  t )  II  M  =  sup  (1+|x|  )“[q(x,t)| 


L  (M)  M<x<°» 

a 


<  sup  ( i+|x| )a|h(x,t) j  +  2  sup  (1+|x|)a  •  /  |h(x+y,t) | |b(x,y, 
M<x<°°  M<X<°°  0 


+  sup  (1  +  |x|)J  /”jh(x+y,t) | | (b*b) | (x,y,t)dy 
M<x<°°  0 


y,t)dy. 


t)  |dy 
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<  lh(*,t)l  ^  +  2  sup  (1+|x|)a  sup  |h(x+y,t)|  •  Y(x,t)eY*X,t* 

L  (M)  M<x<°°  0<y<« 

a 


+  sup  (1  +  | x | ) a  •  sup  |h(x+v,t)|  •  Y2(x,t)e2Y^X,t^ 


M<y<° 


0<y<° 


<  »hf,t)» 


(1+Y(x,t)eY(x,t) )2. 


La(M) 


3cj 

To  estimate  we  note  that 


(5.24)  •|^(x,t)  =  y£(x,t)  +  2  /  |^-(x+v,t)b(x,y,t)dy 


3h, 


+  /"  ^~(x+y ,  t )  (b*b ) (x ,y , t )dy 


+  2  f°  h(x+y,t)-4£-(x,y,t)dy  +  f°  h(x+y,t)  4-(b*b)  (x,y,t)dy 


'It' 


IT 


3h  3b 

The  terms  are  estimated  as  above,  and  the  terms  are  easily  seen  to 

provide  the  remaining  terms  in  the  estimate  (5.21)(iii)  above. 

3q 

Finally,  we  consider  g^-: 


3a  3h  r00  3h 

(5.25)  g^-(x,t)  =  g~(  x ,  t )  +2  J  g^-(x+y,t)b(x,y,t)dy 


3h , 


r”  3h 

+  J  g^-(x+v,t)  (b*h)  (x,v,t)dv 


d  V) 

h(x+y,t)  g^-(x,y,t)dy 


+  2  /  h(x+y,t)(b  *  4^-)  (x,y,t)dy. 

o  ‘■’X 
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ah 

Aqain  the  j — (x,t)  terms  lead  to  estimates  with  a  factor 

(1 

3b 

The  remaininq  terms  (those  involving  j^-)  lead  to  bounds  as  follows: 
(5.26)  /  (1+|x|)a/  |h(x+y,t) | )|^(x,y,t) |dydx 


+/  (1+|x|)“/  |h(x+y,t) | |b  *  |^(x,y,t) |dyd 


<  /°°(1  +  !x|)a/  jh(x+y,t)  |{  |q(x+y,t)  | 

M  0 

+  eY ^X,t^n <x+y,t)n (x,t)}dydx  •  (1  +  sup  0(x,y,t)|) 

y 

<  r  (1+|x|)a  sup  |h(x+y,t) |{  f°  |q(z,t) |dz 

M  0<  y<“  x 


+  eY  ^X,t\ (x,t)b(x,t)}  •  ( 1 +n  (x,  t  )eY  *X' ^  )dx 


/  (1+|x|)a  sup  |h(x+y,t)|  •  n (x,t) (1 +Y  (x,t)eY ^X,t^ ) 
M  0<  yO 


•  (1-wi(x,t)eY(x,t>)dx 


<  "h(-,t)H  oo(1+Y(M,t)eY(M,t))n-Kl(M,t)eY(M,t))  •  Y  (M,  t ) . 

L 

a 

This  completes  the  proof  of  Lemma  5.4. 

We  shall  need  estimates  on  the  difference  between  two  potentials 
q^  1  *  (x  )  ,q(  (x )  qiven  by  (5.3)  from  respectively. 

Lemma  5.5.  Let  b^.b^)  he  as  in  Lemma  5.4  above  and  satisfy  the  estimates 
(5.15)  (i)-(iv).  Then  if  ^  are  given  by  (5.3)  with  b  =  b^^,  we  have: 


1 


(5.27) 


(i) 

llq(1)C,t)  -  q<2)C,t)ll  < 

La(M) 

C  #h(*,t)»  «Y(M,t) 

La(M) 

(ii) 

Nq(1)C,t)  -  q(2)C,t)ll  „(M)  < 

L 

a 

CllhC  , t ) II  „  (M,t) 

La(M) 

(iii  ) 

a„(i>  aj2) 

tll[> 

l\m) 

a 

+  (C  6u (M,t)  •  #h(* ,t)«  + 

La(M) 

C  6n (M,t)  +  C,  6 Y (M  #  t ) ) 

3  4 

(iv) 

a  (1)  a  (2) 

— e,t)  — 2 — (*,t)# 

l’(M) 

3  h 

<  C  W-,t)l  «Y(M,t) 

X  L„(M) 

a 

+  CJhC  , t ) D  m  •  6y  (M , t ) 

VM) 

where 

Cj ,  j  =  1,***,5  depend  only  on 

Y^M.t),  nt(M,t),  U4(M,t), 

<Mf  t )  t  =  1,2. 

Proof.  From  (5.3),  we  have 


(5.28)  q(1)(x,t)  -q(2)(x,t)  =  2/  h ( x+V, t ){ 2(b ( 1  * ( x, y , t)  -b(2)(x,y,t) 

0 


,w<1)  *  w(2)  *  .(2>w 

+(b  *  b  -b  *  b  )(x,y,t))dy 


whereupon 


(5.29)  llq^'c.t)  -  qt2)(*,t)H  1  <  J  (1  +  |x|)a  /  |  h  ( x+y ,  t )  |  {  2 1  b  * 1  b^2) 


L  (M)  M 
a 


+  |b(4)  <  bl4)  -  b(2)  *  b( 2) | } dvdx 


<  2»hC,t) 


M  Y.(x,t) 
r"  1 

J  e 


L  (M)  M 
a 


Y  -(x,t) 

( 1 +Y  ^  ( x ,  t )  e  )6n(x,t) 
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I 

I 


+  »hC  ,t)»  J  (e  n^x.t) 

L  (M)  M 
Ct 


Y-(x,t)  Y  (x,t) 

+  e  2  n  (x,t))e  (1+Y2(x,t)e 


Y  _(x,t) 


\  Af  W  2 

n+Y_(x.t)e  )6n  (x,t)dx 


Y ,  (M,t) 


Y „(M,t) 


<  Hh(*,t)»  e 

L  (M) 
a 


1  *  ( 1 +Y 2 (M , t)e  2  )6y (M,t) [ 2  + 


Y,(M,t)  Y2(M,t) 

(e  1  D (M,t)  +  e  n2(M,t)  )]. 


The  same 


idea  for  lT(M)  leads  to  (ii)  above.  Differentiating  (5.28)  with 


respect  to  t,  we  have: 


(5.30) 


3t  C  I/(M)  1*^  CM) 


+  2  r  (1  +  U|)a  /“  |h(x+y,t)  ( ( Ib^1  b£2)^ 

M  0 


*  |(b"1-  b<2>,  •  (b''’*b‘2’)l  *  K'’-<2’ 


where 


(b(1)  +  b(2))|)dydx  <  b|^-C  ,t)H  •  C  +  2»h(* ,t)#  1  *  D 

v  t  •  *  f  >1  v  T  f  M  ^ 


La(M) 


La(M) 


n  =  fObf^Cx^rt)  -  bJ.z,(x,-,t)B  m ( i  +it b ( 1  >  +  h(2)(x,-,t)H  1f 

t  T  L 


+  llbn)(x,-,t)-b(2)(x,-,t)ll  Jb^’  +  b[2)(x,-,t)»^}dx 


Y  (M,t)  Y_(M,t)  Y  ^  (Mr  1 5 

<  e  (1+Y2(M,t)e  )(HY,(N(t)e 


Y  _(M , t )  Y  (M,t) 

+  Y  (M,  t)e  2  >{«W(M,t)  +  e  Y.(M,t)6n(M,t)} 

2  ^ 


Y  ,(M,t,  ) 


(  1+Y  (M,t)e 


Y  2(M,  t)  )(SY(Mft) 
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l*1  (Mf  t)  ( 1  +  Y1  (M,  t)e 


(e 


Y^M.t) 


Y^M.t) 


Y2(M,t)  Y,(M,t) 

+  e  U  (M,t)( 1+Y2<M,t)e  ). 


Differentiating  (5.28)  with  respect  to  x,  we  have: 


(5.31) 


3q 


(1) 


(x,t)  -  — (X|t)  =2  /  -^-(x+y,t)  [bu  bl^J  (x,y,t)]dy 


(2) 


3h 


(D  _(2). 


+  f°  h(x+y,t){2(b^  *-  b^2) )  +  2(b(1)*  b^U  -  b(2)*  b^2))}dy. 


Estimating  as  in  Lemma  5.4(iv),  we  obtain 


(5.32) 


a  (i)  a,J2)  aiT 

'“fc- (*'t)  1  <C«g7(x,t)« 

l!(M)  X  lNm) 

a  a 


+  2«hC,t)l  w  /*dx  /°“dy(  |b^  )-b^2)  |  +  |b(1)*  b*1*  -  b<2)*  b 
L  (M)  M  0 


X  X 


<2! 

x 


<  c4V,t)l  +  2#h(*,t)l  w  /°dx{»b(1)(x,*,t)  - 

r  W»i\  r  /  V  U  ^ 


L  (M) 
a 


L  (M)  M 
a 


(2) 


+  »(b(1)  +  b(2))(x,*,t)«  <b(1)  -  b(2)(x,*,t)l  > 

X  X  L1  L1 


,3h , 


<  Clg-C,t)«  +  2E»h(n,t)»  M  «Y(M,t) 

X  l\m)  L  (M) 

a  a 


where 
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IT  T  T  T 

B  5  (l  +  Y,e  1  +  Y2e  2 ) { ( 1 +Y 1 e  1 ) ( 1 +Y2e  2 )  +  Y f  2\ 


+  e  '(1+Y2e  2)(Y1(1+Y1)  +Y2(1+Y2)) 


(for  notational  ease  we  have  omitted  the  arguments  of  Y^»Y2  which  are 
M,t).  This  follows  from  a  straightforward  estimate  on  b^1^  -  b^2^  and 
completes  the  proof  of  the  lemma. 

We  are  now  ready  to  prove  Theorem  5.1,  as  we  have  in  fact  derived 
Lipschitz  estimates  for  the  mapping  q  (“►  $(q),  whose  Lipschitz  constants 
depend  linearly  on  various  norms  of  h(x,t). 

Proot  of  Theorem  5. 1 .  Let  J  be  a  compact  subinterval  of  {t:  0  <  t  < 

We  define  the  Banach  spaces  BM  of  functions  via  the  norm 


M 

|q|  I  L  s  sup  l  UaVqc.t,! 

tGJ  s+3r<6  L 


+  H9r3Sq(  •  ,tll 
t  x 


1 3/4-6-( 3r+s )/2 

1 


(M) 


L4-(3r+s)/2(M) 


where  M  will  be  suitably  chosen,  as  discussed  below. 

By  the  estimates  of  the  above  lemmas  and  their  direct  extensions  to  the 

higher  order  derivatives  contained  in  the  norm  above,  the  mapping  q  I-*  $(q) 

maps  Bm  into  itself  tor  every  M,  Moreover,  if  we  consider  two  functions 

q1 ,q2  in  the  unit  ball  of  (which  therefore  bounds  sup  Y(M,t),  n.(M,t), 

tGJ  3  3 

\  Mft  ,  by  1),  then  it  is  clear  that  an  estimate  of  the  form 

j  3 

i  I  ^(q,)  -^q2)lllM<  |  |qn  -  q2l  I  lMl  I  |h|  |  |M 

holds,  where  K  >  1  is  an  absolute  constant  (involving  the  maximum  of 
C,L),E  as  above).  Thus  choosing  M  so  that  11,111 

contraction  and  h  belongs  to  the  unit  ball  in  BM.  Hence  on  the  interval 
[M,t”),  equations:  (5.2),  (5.3)  have  a  unique  solution  q(x,t). 


h|  | |M  <  — ,  $  is  a  strict 
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To  complete  the  proof  of  the  Theorem  5. i ,  we  must  show  that  this  solution 
extends  to  all  finite  x  and  that  as  t  +  0,  q(x,t)  ♦  Q(x)  in  X3_5*  The 
continuation  argument  follows  that  of  {7]  and  proceeds  as  follows: 

Suppose  that  a  solution  of  (5.2),  (5.3)  in  BM  exists  only  up  to  some 
finite  value  of  M  -  i.e.  M*,  the  mfimum  of  ail  M  for  which  a  solution 

v  * 

exists  for  x  >  M  and  belongs  to  BM,  is  finite.  Then  for  M  =  M  +  e,  we 
have  a  solution  q(‘,t)  6  B  for  every  e  >  0.  There  is  a  corresponding 

b(x,y,t)  via  (5.2).  Extend  b(x,y,t)  for  x  >  M  -  E  as  the  constant  (in 

x)  function  b(M,y,t).  Since  h(x,t)  exists  for  all  x,  we  can  consider 

*  * 

(5.2),  (5.3)  on  the  x-interval  (M  -  e,M  +  e).  By  arguing  as  in  the  above 

lemmas,  it  is  easy  to  show  that  our  system  has  a  unique  solution  for 

b(x,y,t)  in  a  ball  relative  to  b(M,y,t)  in  a  space  with  norm: 

§up  A  l  b(x,y,t)»  +  1  »  „  } 

xS(M  — c ,M  +e]  3r+s+o<6  y  L,  L 

t6J  C<1  lyl  (yl 

for  £  sufficiently  small.  The  corresponding  q(x,t)  exists  for 

* 

x  >  M  -  e  and  belongs  to  B  ,  which  contradicts  the  definition  of 


*  * 

M  .  Thus  M 


* 

M  -e 


The  proof  that  q(x,t)  ♦  Q(x)  in  L3  g(M>  follows  for  M  sufficiently 
large  by  iteration  and  for  all  semi-infinite  intervals  by  continuation  as 
above.  This  completes  the  proof  of  Theorem  5.1. 

We  remark  that  analogous  results  hold  for  initial  data  in  with  the 

same  proofs.  Note  also  that,  as  an  immediate  corollary  of  Theorem  5.1,  we 
have  the  result  that  with  L1  initial  data  of  compact  support,  the  solution 


is  C  for  t  >  0  (see  [6]). 
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Verification  of  the  KdV  Equation 


In  order  to  prove  that  our  generalized  solution  constructed  above  is 
indeed  a  solution  of  the  KdV  equation,  it  is  sufficient,  as  remarked  by  Tanaka 
[19],  to  prove  the  following  result 

Theorem  6.1.  bt(x,y,t)  +  bxxx<*,y,t)  -  3q(x,t)bx(x,y,t)  =  0. 

Proof.  The  idea  is  to  show,  using  the  Schrodinger  equation  and  the  trace 

formula,  that  bt  +  bxxx  and  3q  bx  satisfy  the  same  linear  integral 

equation,  which  implies  they  are  equal. 

Differentiating  the  Volterra  equation 

(6.1)  b(x,y,t)  =  /  q(w,t)dw  +  /Y  /  q(w, t)b(w,z,t)dwdz 

x+y  0  x+y-z 


we  obtain: 


(6.2) 


b  (x, y , t )  +  b  (x,y, t) 
t  xxx 


/  (q  (w,t)  +  q  ( w, t ) )dw 

t  WWW 

x+y 


+  /Y  P  [q  ( w,  t )  (bfc  (w,  z ,  t )  +  b^lw.z.t)) 

0  x+y-z 


+  (qt(w,t)  +  qwww(w't)  )b(w,z,t)  +  Sfq^w.tjb^w.z.t)  >w)dwdz. 

The  trace  formula 

(6.3)  q(x, t)  =  h(x,t)  +2  /  h(x+y,t)b(x,y,t)dy 

0 

+  P h(x+y,t)  /Yb(x,y-z, t)b(x,z, t)dzdy 
0  0 


when  differentiated  (recalling  hfc  +  hxxx  =  0)  yields 


t 


(6.4) 


q  (x,t)  +  q  (x,t)  =  2  /  H(x+y, t) (b  (x,y, t)  +b  (x,y,t))dy 

t  XXX  q  C  XXX 

+  2  /*h(x+ y,t)  •  Jyl b(x,y-z,t)(bfc(x,y,t)  +  b (x,y, t) ) 

0  0 

+  3  bx(x,y-z,t)bxx(x,y,t)]dzdy 

+  6  4—  {/  h  (x+y,t)[b  (x,y,t)  +  /yb(x,y-z,t)b  (x,z,t)dz)  •  dy}. 
x  0  0  X 


The  estimates  of  Section  5  ensure  the  convergence  of  all  these  integrals. 
Combining  (6.2),  (6.4),  we  obtain  the  following  integral  equation  for 

*(*,y,t)  =  +  bxxx(x'y,t): 

4»(x,y,t)  =  2  /  [  (h  *  \^)(w,t)  +  (h  *  (ij/  *  b))(w,t) 

(6.5)  X+y 

+  3(h  *  (b  *  b  ))(w,t)  +  3(h  *  (b  +  b  *  b))  (w,t)]dw 

w  ww  w  w  w  w 

+  /Y  r  (q(w,t)iji(w,z,t)  +  3(qw(w,t)bw(w,z,t)  )w 
0  x+y-z 

+  2  b(w,z,  t)  { (h  *  ij»)(w,t)  +  (h  *  (i|»  *  b) )  (w,  t) 

+  3(h  *  (b  *  b  ))(w,t)  +  3(h  *  (b  +  b  *  b  ))  (w,t)}]dwdz. 
w  ww  w  w  w  w 


Here  we  use  *  to  denote  both  types  of  convolutions  in  (6.4)  above.  Thus  it 
suffices  to  show  that  3  q(x,t)  •  bx(x,y,t)  is  also  a  solution  of  (6.5). 
Making  this  substitution  and  dividing  by  3,  we  wish  to  show: 
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(6.6) 


q(x,t)b  (x,y,t)  =  2  /*  C(h  *  qbw)(w,t)  +  (h  *  (qbw  *  b))(w, 


+  (h  *  <bw  *  bww))(w't)  +  (hw  *  <bw  +  b  *  bw} )w(w,t) ]dw 


+  fY  r  [q2(w,t)bw(w,z,t)  +  (qw(w,t)bw(w,z,t) 
o  x+y-z 

+  2b (w, z, t)  { (h  *  qbw)(w,t)  +  (h  *  (qb^  *  b))(w,t) 

+  (h  *  (b  *  b  ))(w,t)  +  (h  *  (b  +  b  *  b  )  )  (w,  t) }  ]dwdz 

W  WW  W  w  W  w 

From  here  on,  the  proof  is  a  matter  of  direct  calculation,  using  the  following 


ingredients: 


(i)  If  $(x,y,t)  6  C  ,  $  +  0  as  x  ++  ",  then 


♦tx.y.t)  =  ♦(x+y.O.t)  +  /y  r  ~  *wz)(w'z't>dwdz 


0  x+y-z 


(6.7) 


(li)  b  (w,z,t)  =  q(w,t)b(w,z,t)  +b  (w,z,t) 

WW 


and  bw<w,0,t)  =  bz(w,0,t)  =  -q(w,t) 


( iii )  /yu(z)v  (y-z)dz  =  -u(y)v(0)  +  fY  u(z)v(y-z)dz}. 

w  aY  n 


Using  (6.7) (i)  to  re-express  the  left-hand  side  of  (6.6)  as 


(6.8)  -q2(x+y,t)  +  /y  r  I(q(w,t)b(w,z,t))ww  -  (q(w, t )b(w, *, t ) >wz )dwdz 

0  x+y-z 

and  using  (6. 7)(ii)  to  eliminate  terms,  (6.6)  is  equivalent  to: 

(6.9)  Kqq, )  =  Mp) 

where,  for  any  function  g(w,t),  we  define 

(6.10)  I (g (w, t ) )  i  1“  g(w,t)dw  +  JY  /“  g ( w, t )b ( w, z , t )dwdz , 

x+y  0  x+y-z 


and  where  p(w, t)  is  given  by 
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(6.11) 


b)  (w,t) 


p(w,t)  H  q(w,t)[(h  *  b  )(w,t)  +  (h  * 

w  w 


+  2(h  *  b  *  b  )(w,t)  +  (h  *  b  *  b)(w,t)] 
w  w 


+  (h  *  b)(w,t)  +  (h  *  (b  )(w,t) 
ww  w  wz 


+  (h  *  (b  *  b  ))(w,t)  +  (h  *  (b  *  b))(w,t) 
ww  w  w  wz 


+  (h  *  b  *  b  )(w,t)  +  (h  *  b  *  b  )(w,t). 
WWW  w  wz 


Clearly  (6.9)  holds  if  qq^  =  P.  Now,  the  terms  in  p  which  do  not  contain 
the  factor  q  may  be  grouped  pairwise.  For  the  first  pair,  we  have: 

(6.12)  f°  [ h  (w+y , t ) b  ( w , y , t )  +  h  (w+y,t)b  (w,y,t)]dy 

ww  w  w  wv 

0 


=  h  (w+y,t)bw(w,y,t)  =  hw(w,t)q(w,t). 

1  Y=0 

The  second  pair  is: 

j  Ih  (w+y,t)  •  /y  b(w,z,t)b  (w,y-z, t)dz 


(6.13) 


+  h  (w+y, t)  /y  b(w, z, t)b  (w,y-z,t)dz)dy 

w  0  wy 


/  (h  (w+y, t)  /y  b(w,z,t)b  (w,y-z,t)dz)  dy 
w  w  y 

0  0 


+  /  h  (w+y, t)b(w,y,t)q(w,t)dy 


=  q(w,t)  f°  h  (w+y,t)b(w,y,t)dy  by  (6.7)(iii). 


A  similar  reduction  on  the  last  pair  of  terms  in  p  implies 

(6.14)  P(w,t)  =  q ( w, t )  [  h  (w,t)  +  2 ( h  *  b  )(w,t)  t  2 ( h  *  b)(w,t) 

w  w  w 


+  2(h  *  b  *  b  )(w,t)  +  (h  *  (b  *  b))(w,t)] 
w  w 
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which  by  the  trace  formula  (6.3)  implies  p  =  q  q  .  Thus  b„  +  b  -  3  qb 

w  t  XXX  X 

=  0,  and  our  generalized  solution  of  the  KdV  equation  is  indeed  a  bona  fide 
solution. 
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